Quantum critical scaling in magnetic field near the Dirac point in graphene 
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Motivated by the recent measurement of the activation energy at the quantum Hall state at the 
filling factor / = 1 in graphene we discuss the scaling of the interaction-induced gaps in vicinity of 
the Dirac point with the magnetic field. The gap at / = 1 is shown to be bounded from above by 
E(l)/C, where E(n) = vf\/2uB is the Landau level energy and C = 5.985 + 0(1/N) is a universal 
number. The universal scaling functions are computed exactly for a large number of Dirac fermions 
N. We find a sublinear dependence of the gap at the laboratory fields of 10T < B < 50T for realistic 
values of short-range repulsion between electrons, and in quantitative agreement with observation. 



I. INTRODUCTION 

When placed in a magnetic field, graphene exhibits a 
series of incompressible quantum Hall states at filing fac- 
tors / = ±(4n + 2), which are a direct consequence of 
the Dirac nature of its' quasiparticles [l|, 0. Whereas 
this main sequence of states can be understood within a 
picture of essentially non-interacting electrons [H , [H , [j| 
the additional incompressible states at other even filling 
factors, / = 0, and at / = ±1, which become discernible 
at higher magnetic fields [6|, call for additional consid- 
erations. The observation of the incompressible state at 
/ = ±1, in particular, implies a complete removal of 
the fourfold spin and valley degeneracy near the Dirac 
point, and has been hypothesized to be due to electron- 
electron interactions. Several theories of these additional 
incompressible states in graphene have recently been put 
forward @, S, S M HO, 111 • 

The idea of an interaction-induced gap at / = 1 has 
very recently been give n strong support by the experi- 
ment of Jiang et al. [13| which showed the activation gap 
in longitudinal resistivity to be independent of the com- 
ponent of the magnetic field parallel to the graphene's 
plane. Furthermore, gap's dependence on the orthogo- 
nal component of the field appears to be sublinear. Both 
features stand in stark contrast to the behavior at / = 4, 
for example, where the gap scales linearly, and with the 
total magnetic field, which makes it very likely to be due 
to the Zeeman effect Q , [k| . The sublinear field depen- 
dence of the gap at / = 1 in particular conforms to the 
expectation that the interaction-induced gap should scale 
with the Coulomb energy scale, e 2 /el, where / = 1/yB is 
the magnetic length. This natural interpretation, how- 
ever, runs into the following difficulty: with the com- 
monly assumed dielectric constant e ~ 5 the measured 
gap at / = 1 (of ~ 10QK) is more than an order of 
magnitude smaller than the theoretical prediction, which 
would require it to be of a similar size as the energy of 
the first Landau level. The observed energy gap at / = 1 
provides therefore an unanticipated intermediate energy 
scale, in between the Landau level separation (~ 1000A') 
and the much lower Zeeman energy (~ 10K). The origin 
of such an energy scale in graphene in the magnetic field 
is presently unknown. The purpose of the present paper 
is to draw attention to this puzzle and propose a solution. 



FIG. 1: The gap m (in units of vfA) at filling factor 
/ = 1 as a function of the magnetic field B (in units of 
Bo = A 2 ) below the zero- field critical value of the short- 
range coupling g, for N = oo. The top curve corresponds 
to the critical point 5 = (Eq. (1)), and the remaining 
ones to 8 = 0.03,0.07,0.14,0.31,0.7 (top to bottom), with 
8 = (g c — g)/gg c A. The best fit to the experimental data [l3l ] 
(dots) assuming 1/A = 2.5A (i. e. B = 10 4 T) is for 8 = 0.31. 



A possible reason for the smallness of the energy gap at 
/ = 1 may be almost trivial: assuming an order of mag- 
nitude larger dielectric constant, which may be due to 
an accumulated layer of water for example [l4j |. would 
obviously bring the theory and the observation closer 
together. Here we wish to put forward an alternative 
and more general explanation which relies only on the 
short-range effects of the electron-electron repulsion, and 
as such completely avoids the ambiguities in the size of 
graphene's effective dielectric constant. The gist of our 
theory is the following. A purely short-range repulsion, if 
above a critical value, would open a gap in graphene even 
at zero magnetic field [l5| . Such a metal-insulator quan- 
tum phase transition is believed to be continuous. Right 
at the metal-insulator quantum critical point then the 
system becomes scale invariant which, as will be shown, 
implies that upon the introduction of the magnetic field 
the gap at / = 1 behaves as 
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(1) 



vf is the Fermi velocity at the Dirac point, \/2B the 
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energy of the first Landau level in units U = e/c = 1, and 
C is a universal number. We determine the universal 
number in Eq. (1) to be 

C = 5.985 + 0(1/N), (2) 

with N as the number of Dirac fermions. Since 1/N cor- 
rections to simpler universal quantities in the problem, 
such as the critical exponents at zero magnetic field, are 
known to be only a few percent even for N as low as 
two [HI, [IH we expect the above number to be quite 
indicative of the real value of C . At a subcritical value 
of the short-range interactions, likely to be realized in 
real graphene, the sublinear dependence m ~ \f~B crosses 
over to the linear dependence m ~ B at lower fields, with 
the values of the gap at all magnetic fields bounded from 
above by the critical curve in Eq. (1). We compute the 
universal scaling function for the gap dependence on the 
magnetic field and the deviation from the criticality in 
the large-N limit. The family of resulting curves is de- 
picted in Fig. 1. The question posed by the experiment 
is then whether there exists a line with enough curvature 
within the range of the laboratory fields which also agrees 
with the observed magnitude of the gaps. In Fig. 1 we 
argue that the answer is yes. Our conclusion is strength- 
ened further by the inclusion of the long-range part of the 
Coulomb interaction, which provides the leading (log- 
arithmic) corrections to scaling at the large-N critical 
point. We find its effect to be an additional curvature 
to the universal function, and its inclusion thus always 
improves the comparison with the available experimental 
data ( see Fig. 3). 

A necessary condition for the applicability of our the- 
ory is that the external magnetic field is low compared to 
the characteristic lattice magnetic field scale, or equiva- 
lently, that the magnetic length is much longer than the 
lattice constant. Even at the highest fields of ~ 45T 
this condition is comfortably satisfied. Provided that the 
system is not too far in the coupling space from its' quan- 
tum critical point, in the next section we show that in 
general under this condition a single relevant coupling 
constant needs to be included, and a simple scaling law 
in the magnetic field emerges. In the section III we ex- 
plicitly compute the scaling function in the large-N limit 
of the theory, and show how in spite of the complexity of 
the self-consistent equation that determines the gap, this 
is possible to do essentially analytically. In this section 
we also show how the ultraviolet (UV) cutoff drops out 
of the universal scaling function for the long magnetic 
length. In sec. IV leading corrections to scaling due to 
the long-range tail of the Coulomb interaction are in- 
cluded. Finally, we discuss various aspects of our results 
in the concluding sections. 

II. SCALING 

Let us begin providing the details behind the above 
results. To be specific, we assume the spin degeneracy 



to be removed by the Zeeman splitting, and that the 
chemical potential is close enough to the Zeeman-shiftcd 
Dirac point so that / = 1. Assuming further a sim- 
ple quantum critical point at B = we retain a sin- 
gle relevant short-range coupling constant g. We will 
also include the coupling representing the ~ 1/r tail of 
the Coulomb interaction, which remains unscreened in 
graphene: A = 2ire 2 /ev F [ljfl. At "weak" magnetic fields, 
at which l/a ^> 1, one is at liberty to use the contin- 
uum field-theoretic description, with the underlying lat- 
tice entering only at energies above the ultraviolet cutoff 
A = l/a. The internal consistency of such a description 
requires that if another value of the cutoff, say A/6, is 
chosen, the gap in the spectrum m satisfies 

m/A = b- 1 v F (b)F±(\6(b)\,\(b),b 2 B), (3) 

where the functions v F (b), 5(b) and X(b) are such that 
the measurable value of m is independent of the arbi- 
trary factor b. Here we defined a dimensionless parame- 
ter S = (gA)^ 1 — (g c A) _1 for later convenience, with g c as 
the critical value of the short-range coupling. The (engi- 
neering) scaling of the magnetic field in the last equation 
follows from gauge-invariance. The two functions F+ and 
F_ refer to S > (g < g c ) and S < (g > g c ), respec- 
tively. 

We may then choose b = (Bq/B) 1 / 2 = l/a, with B$ = 
1/a 2 , and write 

m = r 1 v F {l/a)F±[6(l/a),\(l/a),B }. (4) 

If l/a ^S> 1 the flow of the couplings is essentially deter- 
mined as at zero magnetic field. Omitting all but the 
single relevant parameter S from the outset becomes jus- 
tified, as all the other couplings are much smaller at a 
large (magnetic) length scale, provided that in the cou- 
pling space the system was not too far from the critical 
point. Assuming a simple critical point at B = near 
which \S(b)\ = |<5|6 1 ^ ; x(b) = A*, and v F (b) = v F b x -\ 
where v and z are the usual correlation length and dy- 
namical critical exponents, we may finally write 

m/(v F A) = (a/l) z G ± [l/^X% (5) 

with £ = a|<5| _iy , as the correlation length. This is the 
universal scaling form of the gap at low magnetic fields, 
dictated by the critical point at S = 0, A = A*, and B = 
[III. The last expression is analogous to the finite-size 
scaling, with the magnetic length playing the role of the 
system's size [l9j . 

Before proceeding with the computation of the univer- 
sal scaling function, let us consider its' limits first. At 
the critical point, 

G_[0,A*] = G+[0,A*]. (6) 

In particular, for z = 1 [l5| this yields the announced Eq. 
(1). At x 1, we must have 

G + {x,\*]~x z - 2 , (7) 
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so that for g < g c one finds m ~ B, as appropriate to 
"magnetic catalysis" (2(j. For g > g c , on the other hand, 
at a; ^> 1, 



G_[x,A*] 



(8) 



which yields m ~ |<5| Z1/ above the critical coupling and 
at B — 0. This is the familiar scaling of the gap near a 
quantum critical point [2lj . 



III. LARGE-N CALCULATION 



+ 



Let us now take A = 0, and compute the scaling 
function G+(x,0) in the limit of large number of Dirac 
fcrmions. The discussion of the effects of long-range 
interaction will be presented in the next section. The 
ground state energy for N species of four-component 
Dirac fermions in magnetic field and when N — > oo was 
derived previously in [10( : 



E(m) - E(0) _ m 2 
N ~ 4^7 



B 



ds 



47T 3 / 2 J S 

[p + 2K(sA 2 )(coth(sB) - 1 



1) (9) 



m is the gap to be determined by the minimization of 
E(m), and g is the dominant among the quartic couplings 
that represent the short-range part of Coulomb repulsion 
flU ]. K(x) is the cutoff function introduced to sum over 
n 7^ Landau levels, which satisfies K(x — > 00) = 1 and 
K(x — > 0) = 0, but is otherwise arbitrary. The param- 
eter p = 2 for / = 0, when the zeroth Landau level of 
each Dirac fermion contributes to the energy difference 
E(m) - E(0), and p = 1 for / = N/2, when half of 
the zeroth Landau levels do not [lOj]. E(m) may also be 
understood as the variational ground state energy of elec- 
trons with an on-site, or nearest-neighbor, repulsion on a 
honeycomb lattice and in low magnetic field, in Hartree 
approximation. 

Hereafter we set p = 1, which would correspond to 
the filling factor / = 1 for the physical case of N = 2. 
Minimizing E(m), after some transformations the gap 
equation can be written in a compact form as 



y = f(y) + 25(yB /B) 1 / 2 , 



(10) 



where the variable y — B/m 2 , and S is as defined below 
Eq. (3), with 



1 

9c 



A 



The function f(y) is defined as 



2 f 00 dt v( tyB 



2yte' 



V 



(11) 



(12) 



The ultraviolet cutoff enters the gap equation by provid- 
ing the scale for the magnetic field (Bo). In the limit 



FIG. 2: The universal scaling function G+(x, 0) at iV = 00. 



of weak magnetic field, B/Bq <c 1, we may replace the 
cutoff function K(x) in the last equation with unity. In 
this (continuum) limit the dimensionless variable y be- 
comes a universal function of the ratio 5(Bo/ ' B) 1 / 2 . The 
only remaining dependence on the cutoff then is in the 
value of the critical point g c which is therefore, as usual, 
non-universal. 

The gap equation may be solved essentially analyti- 
cally for S > by noticing first that at S = the (numer- 
ical) solution is at yo = 17.913. This immediately yields 
the result in Eqs. (1) and (2), with C — \/2yo~. Since at 
S > the solution will be at y > yo, one in fact needs 
the function f(y) only for large arguments, where it can 
be expanded as 



f(y)=uy 1 / 2 +vy- 1 / 2 + 0(y-V 2 ), 



(13) 



with u = 4.13031 and v = 1.84723. The difference be- 
tween f(y) and the first two terms in the expansion on the 
right-hand side is less than a percent already for y > 3, 
and decreases further with y. Comparing with the gen- 
eral form in Eq. (5) we may then write the Eq. (10) in 
terms of the universal function G+(x,y) as 



- 1 + (u + 2x)G+(x, 0) + vG 3 + (x, 0) = 0, 



(14) 



with the higher order terms in G+ entirely negligible. 
This way we find G+(0,0) = 0.236, and G+(a;,0) = l/2x 
for x 3> 1, and uniformly decreasing in between. In sum, 
at N — 00 and A = the Eq. (5) becomes 



L z a 



(15) 



with z = 1 and v = 1 + 0(1/N), and the function 
G+(x, 0) depicted in Fig. 2. 

Although it may not be directly relevant to graphene, 
for completeness we also briefly describe some of the re- 
sults when 6 < 0. Since the gap is finite even at zero 
field above the critical coupling, the solution of Eq. (10) 
now lies at y < y . For weak fields, the relevant regime 
is y <C 1, where f(y) = 4 + 2y + 0(y 2 ). This implies, 
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for example, G-(x, 0) — > x/2 when x 3> 1, in agreement 
with the Eq. (8). We thus find m = |£|/2 + 0(<5 2 ) in the 
zero-field limit, and at N = oo. 



IV. LEADING CORRECTIONS TO SCALING 

Let us now turn on a weak long-range interaction A <C 
1. It represents a marginally irrelevant perturbation at 
the large-N critical point 5 = 0, A* = 15], and thus 
provides logarithmic corrections to the scaling law in Eq. 
(5). At b ^> 1, the Coulomb interaction scales as 



A(6) 



8tt 



(16) 



The effect of the long-range Coulomb interaction is to 
break the pseudo-relativistic invariance of the problem 
and renormalize the velocity as [22j , (23j 



v F (b) 



0(A 2 ))ln6). 



The Eq. (4) then yields, for l/a > 1 

m= (^ + °( A 2 )W^f^G+(— ,0). 



(17) 



(18) 



Interestingly, the inclusion of such corrections improves 
the agreement with experiment. As an illustration, in 
Fig. 3 we display the fit to 

m = VF VB(l + — In— )G+( — ,0). (19) 
l07T B a 

1 / \f~B* defines the length scale appropriate to the mea- 
sured value of vf = 10 6 m/s, which we will treat as a 
fitting parameter. The best fit for e = 1, for example, 
is obtained for B* = 29T. More data seems needed, 
however, to distinguish between this and the alternative 
forms such as a simple m ~ \/~B considered in ref. 13. We 
hope our theory will stimulate more experimental work 
in this direction. 



V. DISCUSSION 

We have argued that taking into account only the 
short-range parts of the repulsive electron-electron inter- 
action in graphene suffices to provide a qualitative, and 
even semi-quantitative understanding of the observed 
magnetic field dependence of the gap at / = 1. Our main 
assumption is that the value of the dominant short-range 
coupling is below but not too far from its critical value, 
so that the aforementioned magnetic field dependence re- 
sembles itself as right at the criticality at the laboratory 
fields. 

Note that the precise nature of the short-range cou- 
pling, i.e. whether it is on-site or nearest-neighbor in- 
teraction, for example, does not matter for the form of 



B 



FIG. 3: The same as in Fig. 1 at S — 0.3, without (dashed 
line) and with (A/167T = 0.27, (e = 1), thick line) the loga- 
rithmic corrections to scaling in Eq. (19). 



the scaling function in the large-N limit we considered. 
It is only the resulting order parameter [To| that will de- 
pend on this. Any order parameter which breaks the 
chiral symmetry of the Dirac Hamiltonian and lifts the 
valley degeneracy in the magnetic field will lead to the 
same large-N scaling function as obtained here. 1 /TV- 
corrections may depend on the type of the order param- 
eter, however. 

The best fit to the experimental data is provided by 
S = 0.31 (Fig. 1), which places the system relatively far 
from the critical point. On the other hand, assuming the 
simplest Hubbard model for short-range interactions, for 
example, would place the critical point at U c / 1 = 4 — 5 
[3. With t = 2.5eV, the usual estimate Q U » 10W 
is then in a reasonable agreement with the value of d ob- 
tained from the fit. It is in fact not uncommon that a 
simple single-parameter scaling works well even reason- 
ably far from the critical point. A celebrated example 
is provided by the scaling plot for classical liquids [2lj . 
where very good scaling of the data is found even at 
T/T c re 0.5. 

For weak Zeeman splitting the gap at / = will also 
obey the critical scaling in magnetic field, and the critical 
function can be found in the large-N limit by choosing 
p = 2 in the Eq. (9). At / = 0, however, the Zeeman 
gap is in competition with the "mass" -gap studied here, 
and there are reasons to believe that the gap at / = 
in reality may be a pure (albeit interaction enhanced) 
Zeeman gap For this reason we have not discussed 
here the scaling at / = in further detail. 

A test of relevancy of the presented theory would be a 
measurement of the field-dependent gap at the filling fac- 
tor / = 1 on a suspended graphene, for example, where 
the dielectric constant e ?» 1. If the gaps are still of the 
same size as in ref. 13 this would eliminate the alterna- 
tive explanation we mentioned in the introduction, that 
the smallness of the gap is caused by the extra screen- 
ing by the surrounding medium. Alternatively, one may 
want to enhance screening by bringing a metallic plate 
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to the graphene layer. If the gap is indeed caused pri- 
marily by the short-range electron-electron interaction as 
discussed here, this should not alter its size much either. 

VI. CONCLUSION 

We have discussed the general magnetic-field scaling 
of the gap at filling factor / = 1 for electrons interacting 
via short-range repulsive interactions in graphene. Scal- 
ing functions and the critical exponents are given explic- 
itly in the large-N limit of the theory, and the leading 



corrections to scaling deriving from the long-range tail 
of the Coulomb interaction are determined. Our theory 
reproduces the recent experimental data for reasonable 
values of short-range interactions. 
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